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Abstract 



By using some basic calculus of multiple integration, we provide an alternative 
expression of the integral 



in which the minimum and the maximum are replaced with two single variables. 
We demonstrate the usefulness of that expression in the computation of orness and 
andness average values of certain aggregation functions. By generalizing our result to 
Riemann-Stieltjes integrals, we also provide a method for the calculation of certain 
expected values and distribution functions. 
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1 Introduction 

Let a, 6 G M U {— oo, +00}, with a < b, and consider an integral over ]a, 6[" 
whose integrand displays an explicit dependence on the minimum and / or the 
maximum of the variables, that is, an integral of the form 



In this note we provide an alternative expression of this integral, in which the 
minimum and the maximum are replaced with two single variables. When the 
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integral is tractable, that alternative expression generally makes the integral 
much easier to evaluate. For instance, when the integrand depends only on the 
minimum and the maximum of the variables, we obtain the following identity 



p rb rv 

/ /(minxj, max Xj) (ix = n(n — 1) dv f{u,v){v — u)^~'^ du, (2) 

J]aM" Ja Ja 



f (minxj, max Xj) (ix = n(n — 1) / 
and hence, for certain functions /, the integral becomes very easy to evaluate. 



The alternative expression we present for integral ([T]) is given in the next sec- 
tion (see Theorem [3l). The method we employ to obtain that expression merely 
consists in dividing the domain ]a, 6[" into n polyhedra chosen in such a way 
that the minimum and maximum functions simply become single variables. 

This method can be very efficient in the evaluation of many integrals that 
would normally require difficult and tedious computations. As an example, 
consider the variance of a sample x G [a, 6]" from a given population, namely 

1 " / 1 " \ 2 



The average value over [a, 6]" of the variance-to-range ratio function can be 
easily calculated by using our method. We merely obtain 

-/ dx= 6- a . 3 

(o — a)" J[a,b]" maxxj — mmxj 12n 



This note is set out as follows. In Section 2 we state and prove the main re- 
sult. In Section 3 we provide an application of our result to internal functions, 
also called Cauchy means, which can be classified according to their location 
within the range of the variables. A similar application to conjunctive and 
disjunctive functions is also investigated. In Section 4 we show how the di- 
rect generalization of our result to Riemann-Stieltjes integrals enables us to 
consider the evaluation of certain expected values from various distributions. 

We will use the following notation throughout. For any n-tuple x, we denote 
by (x I Xj = u) the n-tuple whose ith. coordinate is m if i = j, and Xi otherwise. 
Also, for any integer n ^ 1, we set [n] :={!,..., n}. 



2 Main result 



In this section we present our main result which consists of an alternative 
expression of integral ([T]). We start with a preliminary lemma, which concerns 
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the particular cases of functions involving either the minimum or the maximum 
of the variables. 



Lemma 1 Let f : ]a, b[^~^^ ^ M. be an integrable function. Then we have 



n ^^i^ 
*eN\{i} 



/ du 

J a 




/ dv 

J a J 


Ma, 



n ^^i- 

ieN\{i} 



Proof. Consider the following ra-dimensional open polyhedra 

Pj := {x G ]a, 6[" : x^ > Xj ^ j} (j G [n]). 

They are pairwise disjoint. Indeed, if x G PjHPk, with j 7^ k, then Xk > Xj and 
Xj > Xk, which is a contradiction. Moreover, the union of their set closures 
covers ]a,6[". Indeed, for any x G]a, there is always j G [n] such that 
Xi ^ Xj for all i ^ j. 

Therefore, for any integrable function / : ]a, b\^^^ M, we have 



/ /(x, min Xj) (ix = / /(x, minxj) cix 



which proves the first formula. The second formula can be established similarly 
by considering the polyhedra 

Qj := {x G ]a, 6[" : Xi < Xj Wi ^ j} (j G M). □ 



Lemma [T] is interesting in its own right since it provides special cases of the 
main result. For instance, by applying the first formula, we immediately obtain 
the following identity, which will be used in the next section (see Example [8]). 
For any C [n], we have 

/ f ( min Xi) dx = {b-a)''~\^\\S\ f{u){b-uy^\~Uu. (4) 

J]a,b[^ ^ «e5 / Ja 

Remark 2 We note that, for bounded and continuous functions f , LemmaUl 
can also be derived from the classical Crofton formula, well known in integral 
geometry (see for instance f^). In the appendix we present an alternative proof 
of Lemma U\ constructed from Crofton formula. 
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Let us now state our main result, which follows immediately from two appli- 
cations of Lemma [TJ 

Theorem 3 Let n ^ 2 and let f : ]a, b\^~^'^ — » M 6e an integrahle function. 
Then we have 



U, Xk = v) Y\_ dxi. 

i<^ln]\{j,k} 

A direct use of this result leads to formula Indeed, as the integrand is 
symmetric in its variables, we simply need to consider 

/(x, U, v\Xj = U, Xk = v) = . . . , Xn-2, v) , 

V — U 

where, for any fixed a < u < v < b, the right-hand side is a quadratic polyno- 
mial in Xi, . . . , Xn-2- 



3 Application to aggregation function theory 

We now apply our main result to the computation of orness and andness 
average values of internal functions and to the computation of idempotency 
average values of conjunctive and disjunctive functions. 

3. 1 Internal functions 

We recall the concept of internal functions, which was introduced in the theory 
of means and aggregation functions. 

Definition 4 A function F : ]a, 6[" —>■ is said to be internal if 
minxj ^ F{x.) ^ maxxj (x G ]a, 

Internality is a property introduced by Cauchy [1] who considered in 1821 the 
mean of n independent variables function F{xi, . . . , Xn) which 

should be internal to the set of Xi values. Internal functions, also called Cauchy 
means, are very often encountered in the literature on aggregation functions. 
Most of the classical means, such as the arithmetic mean, the geometric mean. 



/ /(x, min Xj, maxxj) (ix 

J]a,bl" 

J2 dv du /( 



j,k=l 
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and their weighted versions, are Cauchy means. For overviews on means and 
aggregation functions, see the monograph [2| and the edited book [3]. 



It is straightforward to see that a function F : ]a, 6[" — M is internal if and 
only if there is a function / from ]a, 6[" \ diag(]a, to [0, 1] such that 

F(x) = minxj + /(x) (maxxj — minxj), 

where diag(]a, := {(x, . . . , x) G ]a, 6p : x G ]a, h\}. 

Starting from this observation, Dujmovic [5] (see also [7]) introduced the fol- 
lowing concepts of local orness and andness functions, rediscovered indepen- 
dently by Fernandez Salido and Murakami [9] as orness and andness distribu- 
tion functions. 

Definition 5 The orness distribution function (resp. andness distribution 
function^ associated with an internal function F : ]ci, 6[" ^ R a function 
odfp (resp. adfp), from ]a, 6[" \ diag(]a, to [0, 1], defined as 

F(x)-minxi / n maxxi-F(x) 

odjpl^s.) = (resp. adfpM = )■ 

max X,- — mm x,- max x, — mm x. 



Thus defined, the orness distribution function (resp. andness distribution func- 
tion) associated with an internal function F : ]a, 6[" M measures, at each 
x G ]a, the extent to which -F(x) is close to maxxj (resp. minxj), that is, 
the extent to which -F(x) has a disjunctive (resp. conjunctive) or orlike (resp. 
andlike) behavior. 

To measure the average orness or andness quality of an internal function over 
its domain, Dujmovic [5| also introduced the concepts of mean local orness and 
andness, later called orness and andness average values by Fernandez Salido 
and Murakami [9]. 

Definition 6 The orness average value (resp. andness average valuej of an 
internal and integrable function F : ]a, b["' — >■ M defined as 

^p = orf/Hx) rfx (resp. ^p = adf,{.) d^). 



As an immediate property, we note that 

odfri^x.) + adfp(s) = 1, 

which entails odfp + adfp = 1. Thus, as expected, both odfp and adfp render 
the same information and hence we can restrict ourselves to the computation 
of odfp. 
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Even though the computation of odfp remains very difficult in most of the 
cases, Theorem [3] enables us to rewrite this integral in a more practical form, 
namely 

odfp = -z — / dv [ du [ "^^"^ ^ ^ — — - IT dxi 

ip-aY^'^Ja J a J]uA-^ v-u .eNVO'fc} 



The following two examples demonstrate the power of this formula: 

Example 7 Let us calculate the orness average value over [0, 1]" of the geo- 
metric mean 



I In 



i=l 



The case n = 2 is straightforward. Using ^ with f{u,v) = , we obtain 

odfci^) = ln4 - 1. 

Assume now that n ^ 3. As the integrand is a symmetric function, we can 
simply consider 

G*^"^(x I Xj = u,Xk = v) = . . .,Xn-2,u,v) 

and hence, we have 



/ G*-"''(x I Xj = u,Xk = v) Yl dxi 

i&[n]\{j,k} 



Then, using the binomial theorem and observing that -^j^ = iJ^iZoivY' 
obtain 



dv / du 

Jo V — u 

i=0 k=0 



k J ni + {k + l){n + l) 

EE \ (-1)' / x'''+'''+'+^ dx 
i=o k=o\ ^ J 

pi rjfi-(\ rj.n+l\n-2 



lo 1 — 



dx. 
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Finally, 



odfcw = n{n - 1) — — - / dx 



n+lj Jo n-2 

The values ofodfaM for n = 2,3,4,5 are ln4 - I, ^ - %, ^ - if, 
^^^§(25-1175)-^^, respecUvely. 

Example 8 Let us calculate the orness average value over [0, 1]" of a function 
of the form 

SC[n] 

where the set function a : 2t"'l — ^ M fulfills 



SC[n] 



a(0) = and E ^i^) = 1 

SC[n] 

and is chosen so that the function C^") is nondecreasing in each variable. Such 
a function is known in aggregation function theory as a Lovasz extension or 
a discrete Choquet integral (see for instance fTUSJ^). As particular cases, we 
can consider any weighted mean J2i WiXi and any convex combination J2i WiX(_i) 
of order statistics. 

The case n = 2 is easy. We simply obtain odf ^(2) = ^(^a{{l}) + a({2})^. 

Assume now that n ^ 3. For any ^ u < v ^ 1 and any x e [u, we 
have 

C^"''(x I Xj = = f) = E cl{S)u + E (^{S) a;j + E ^^i^) niinxj. 

S3k S^k 

Setting s := \S\, from @) it follows that 



/ Cj^"'\'x I Xj = u,Xk = v) TT dx. 



i&[n]\{j,k} 
SBk S^k 



\SBj 



Then, we obtain 
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V — U J]u,v[ 



i6[n]\{i,fc} 



Summing over j, k = 1 . . . ,n, with j ^ k, and then rearranging the terms we 
finally obtain 



-Tu / fQ\ n{n-l) + s-l \ 1 

^ (n-l)(« -2)(. + l)j - ^731 



5C[n] 



n(n — 1) + s — 1 1 



{n - l){n - 2){s + 1) 



T? — 1 ^ S + 1 

which includes the case n = 2. 

To overcome the difficulty of calculating intractable orness average values, Du- 
jmovic [6] introduced the next concept of global orness and andness measures 
(see also [7f9] ). Denote by F the average value of any internal and integrable 
function F : ]a, 6[" — >■ M over its domain, that is. 



{b - a)" J]a,b[ 



Definition 9 The global orness value (resp. global andness value^ of an in- 
ternal and integrable function F : ]a, 6[" ^ M zs defined as 



F-Min , , Max-F , 

ornessi? = = — = (resp. andness^ = = — ^=/, 
Max-Min ' Max - Min^ 

where Min and Max are, respectively, the minimum and maximum functions 
defined in ]a, fo^. 

For example, considering the geometric mean G'*^"''(x) = IliLi^i'^" [O5 1]"; 
we simply obtain 

1 n + l—-j 1 + 1 / n 

ornessQ{n) = 1 G*-"-' = h 



n — 1 n — 1 n — 1 n — l\n + 1 

Considering the discrete Choquet integral C^"^ in [0, 1]", as defined in Exam- 
ple [HI we get 
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n — 



\s\ + i 



n + l 



SC[n] 



n — \S 



n — 



5C[n] 



Surprisingly enough, in [0, 1]" we have 



orness„{n) = odf („), 



that is, for any discrete Choquet integral, the global orness value identifies 
with the orness average value, a result already reached by Fernandez Salido 
and Murakami [9] for the special case of symmetric Choquet integrals, that is, 
convex combinations of order statistics. 

The interesting question of determining those internal functions F : — ^ 
M fulfilling the equation ornessi;' = odfp remains open. 

3.2 Conjunctive and disjunctive functions 

Let us now consider conjunctive and disjunctive functions. 

Definition 10 A function F : ]a,6[" ^ W is said to be conjunctive (resp. 
disjunct ivej if 



Prominent examples of conjunctive (resp. disjunctive) functions in the liter- 
ature are t-norms (resp. t-conorms), which are symmetric, associative, and 
nondecreasing functions, from [0,1]^ to [0,1], with (resp. 1) as the neutral 
element. For an account on t-norms and t-conorms, see for instance the book 
by Alsina et al. pj. 

Clearly, a function F : ]a, b[^ — > M is conjunctive (resp. disjunctive) if and 
only if there is a function / : ]a, 6[" [0, 1] such that 



Just as for the orness and andness distribution functions, we can naturally 
define the concept of idempotency distribution function associated with a 
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conjunctive (resp. disjunctive) function F : ]a, — > M as a measure, at each 
X G ]a, of the extent to which F is idempotent (i.e., such that . . . , x) = 
x), that is, the extent to which F is close to minxj (resp. maxx,). 

Definition 11 The idempotency distribution function associated with a con- 
junctive (resp. disjunctive) function F : ^ M a function idfp : 
]a, — > [0, 1], defined as 

F(x) - a h - F(x) 

idfpiy^) = —. (resp. tdfri^) = 7 )■ 

mm Xj — a o — max x,- 



We can now introduce the concept of idempotency average value as follows. 

Definition 12 The idempotency average value of a conjunctive or disjunctive 
function F : ]a,b[^ ^ W is defined as 



[b — a)" J]a,b[" 



According to Lemma [T], for any conjunctive function F : ]a, 6[" —>■ M for 
instance, we can write 

idfp = — — [ du [ ^^'^ ^ — — - TT dxi. 

[b — ar ~{Ja J\uM^-^ u — a -^r 



The following concept of global idempotency value was introduced by Kolesarova 
[11] for t-norms as an idempotency measure: 



Definition 13 The global idempotency value of a conjunctive (resp. disjunc- 
tive) function F : ]a, 6[" ^ M is defined by 

ide,.p, ^ ^ (res,. ide,np, ^ 

Mm — a b — Max 

Example 14 Let us calculate the idempotency average value and the global 
idempotency value over [0, 1]" of the product 



p(")(x) = n 



X. 



1} 



i=l 



which is a conjunctive function. 
We immediately obtain 
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^. n-1 

idf p{n) = n 



du i T\ Xi] dxi ■ ■ ■ dxn-i 

Jo V / 

^ Jo 



2' 

Setting u = f and then using the classical beta function 



B{a,b) = - tf-^ dt, 

Jo 



we obtain 



= iiB(l/2.„) = ^r<l/2'r<"' 

Or), V / ? / Or). 



2" ' ' ' ' 2" r(n + 1/2) 

in— 1 



+ 1 

X = . 

2" 



On the other hand, we have 

„ n 

idempp(„) = (n + 1) / TT (i: 
4 Application to probability theory 



Since the idea behind our results merely consists in breaking the integration 
domain into smaller regions, Lemma[T]and TheoremOcan be straightforwardly 
extended to Riemann-Stieltjes integrals, thus making it possible to consider 
average values from various probability distributions. 

Consider a measurable function g : M"'^^ — *• M and n independent random 
variables Xi, . . . {i G [n]) having distribution function Fi{x). Define 

the random variable Yg as 

Yg := (^(X, minXj, maxXj), 

where X denotes the vector (Xi, . . . , 

The direct generalization of Theorem [3] to Riemann-Stieltjes integrals can be 
used to evaluate the expected value of Yg, namely 

^[^g] = / fi'(x, minxi,maxxj) dFi(a;i) ■ ■ ■ciF„(x„). 
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It can also be used in the evaluation of the distribution function of Yg, which 
is defined as 

Fgiz)=mi^-y,)] 

= / h(z - g{:s.,mmxi,maxxi)) dFi{xi) ■ ■ ■ dFn{xn), 

where if : M — >• {0, 1} is the Heaviside step function, defined by H{x) = 1 
if X ^ 0, and otherwise. Note that the case where Yg is a lattice polyno- 
mial (max-min combination) of the variables Xi, . . . , X„ has been thoroughly 
investigated by the author in p^lIT^ . 

To keep our exposition simple, let us examine the special case where Yg de- 
pends only on minXi and maxXj, that is, 

Yg ■= 5((minXj,maxXj), 

where : — R is a measurable function. In this case, our method immedi- 
ately leads to 




j,k=l 



In the particular case where the random variables Xi, . . . , are independent 
and identically distributed, each with distribution function F{x), the expected 
value clearly reduces to 

/oo ™ , , n— 2 

dF{v) g{u,v) (F{v) - F{u)) dF{u), (5) 
-oo J—oo ^ ^ 

which generalizes formula ([2]). 

Example 15 For exponential variables Xi, . . . , Xn, each with distribution func- 
tion F{x) = 1 — e^'^^ (x > 0), we simply have 

E[Yg] = nin - 1) J \e-^''dvj ^(m, t;) (e"^" - e"^^) Xe-^^'du. 

Using the change of variables x = e~^'^ and y = e^'^" — e^'^'" , this integral can 
be easily rewritten as 

E[Yg] = n{n - 1) dx £ y^~^ g{ - j ln(x), -j ln(x - y)) dy. 
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Example 16 Let us calculate the distribution function and the raw moments 
of the random variable 

^ max Xi — mill Xj 
maxXj 

from the uniform distribution over ]0, 1]". 

The raw moments can be calculated very easily from For any integer r ^ 0, 
we have 

Eirn = n(n-i) t dv r f^^^y (v-ur-^du= """^ . 

^ ^ ^ 'Jo Jo \ V J ^ ' n + r-1 

On the other hand, the distribution of Y is simply given by 

that is, 



F(z) =n(n-l) I dv I H ( z - ^—-^) {v - uj'-^ du 



F{z) 



0, ifz^O, 

n{n-l) [ dv r {v - uf"^ du = z''~\ ifQi^z^l, 

Jo Jv{l-z) 

n{n-l)( dv r iy - uf-"^ du = I, ifl^z. 
Jo Jo 



Appendix: The use of Crofton formula 

We provide a proof of Lemma [1] as a direct consequence of Crofton formula. 
See [8j for a very good expository note on Crofton formula. 

For < V < V + h < V , let D{v) he a domain of area or volume v. By a domain 
we mean a closed bounded convex set in M*^ for some k. Assume that for vi < V2 
we have D{vi) C D{v2). Let Ai, . . . , A„ be n independent points randomly 
selected with uniform distribution in D{v + h) and let Y = /(Ai, . . . , A„), 
where / is a bounded function. Let A{v) be the event that all the points are in 
D{v) and let Bj{v, h) be the event that Xj e D{v + h) - D{v) and Aj G D{v) 
for all i J. Let = E[Y\A{v)] and let h) = E[Y\Bj{v, h)]. 

In its nonsymmetric version, Crofton formula states that, if lim/j^g h) = 
^j{v) exists and is continuous for all j, then for V > 0, 

E[Y] = ^,{V) = —Y. / v-'^i,{v)dv. 

V JO 

Choosing V = h — a and D{v) = [a,a + v] (which implies D{V) = [a, b]), we 
simply obtain 
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^^^^ = [ , ^ W 

(O — a)" J[a,b]" 



ra+v 


pa+v+h 




/ (ixi ■ ■ 


■ / dXj ■ ■ 


■f 


la 


J a+v 


J a 



^ f f{x\xj = a + v) n dxi. 



«6N\{i} 



If / is continuous in each argument then hm/j^o l^'^i'^^ h) = ^j{v) exists and is 
continuous. According to Crofton formula, we obtain 



I 1^ r-V r 

^liV) = jT^Y. f{^\xj = a + v) n dxi 



which proves the second formula of Lemma [H The first one can be established 
similarly by considering D{v) = [b — v, b]. □ 
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